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ABSTRACT 



Goertz and Morfill (1983) propose that spokes on Saturn's rings form under 
radially moving plasma clouds produced by meteoroid impacts. We demonstrate 
that the speed at which a plasma cloud can move relative to the ring material 
is bounded from above by the difference between the Keplerian and corotation 
velocities. The radial orientation of new spokes requires radial speeds that are at 
least an order of magnitude faster. The model advanced by Goertz and Morfill 
fails this test. 

Subject headings: Planetary rings, Saturn — Saturn, magnetosphere 
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1. Introduction 

The nature of the "spokes" in Saturn's rings remains a matter of speculation 25 years 
after their discovery by the Voyager spacecraft. A brief summary of their properties is given 
here; for more details see Mendis et al. (1984) and references therein. 

1. Spokes are transient radial albedo features superposed on Saturn's rings. 

2. Spokes are composed of dust with a narrow size distribution centered at 0.5 micron. 

3. Spokes have optical depths of about 0.01. 

4. Spokes are only seen near corotation, which is where the Keplerian angular velocity of 
the ring particles matches the planet's rotational angular velocity. Corotation occurs 
in the outer B ring. 

5. Individual spokes measure about 10,000 km in radial length and 2,000 km in azimuthal 
width; they extend over about 10% of the ring radius. 

6. Spokes are seen preferentially on the morning ansa of Saturn's rings, and are most 
closely radial there. 

7. Spokes fade and are distorted by differential rotation as they move from morning 
toward evening ansa. 

8. A few observations have been interpreted as showing the birth of individual spokes 
within 5 minutes along their entire lengths. This timescale implies a propagation 
velocity of at least 20 km s -1 . 

9. Spokes are only observed at small ring opening angle to the Sun (McGhee et al. 2005). 

The above observations suggest that spoke formation involves the sudden lifting of a radial 
lane of dust grains from the surface of the rings. Their subsequent fading and distortion is 
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compatible with the elevated dust grains moving on Keplerian orbits that intersect the ring 
plane half an orbital period (i.e. about 5 hours) later. 

Currently the most popular model for spoke formation is that of Goertz and Morfill 
(1983, GM). GM propose that the formation of a spoke is initiated when a meteoroid 
impacts the ring and creates a dense plasma cloud. Electrons from the cloud are absorbed 
by the ring producing a large electric field which levitates negatively charged dust grains. 
The grains enter the cloud where they absorb additional electrons. Overall charge neutrality 
is maintained by the net positive charge of the plasma. 

Because the dust grains are massive, they move on Keplerian orbits. The plasma in 
which they are immersed is however tied to the magnetic field lines which pass through the 
ionosphere of Saturn. The motion of the negatively charged dust relative to the positively 
charged plasma produces an azimuthal electric field that causes the plasma cloud to drift 
radially. GM argue that the plasma cloud will continue to levitate dust grains as it moves. 
According to their calculations, the plasma cloud drifts in the radial direction, away from 
corotation, at 20 - 70 km/s. This velocity is sufficient to account for the formation of a 
radial spoke of length 10,000 km within 5 minutes. 

We perform a self-consistent calculation of the plasma cloud drift velocity in §2. It 
establishes that the drift velocity cannot exceed the difference between the local Keplerian 
and corotation velocities. This upper limit is of order 1 km s" 1 in the region where spokes 
are observed. We reveal the source of GM's error in §3 and estimate the correct drift 
velocity in §4. In §5 we comment on the response of Morfill & Thomas (2005) to these 
issues. A short summary in §6 concludes our paper. 
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2. Drift of a Plasma Cloud 

The essence of GM's model is illustrated in Fig. 1. The part of the ring plane at the 
base of the plasma cloud has a finite Hall conductivity: in the presence of an electric field 
in the local rest frame of the ring particles, the net positively charged plasma will drift 
relative to the negatively charged dust grains which are embedded within it. Away from 
the plasma cloud the ring plane has negligible electrical conductivity. 

Overall charge neutrality of the dusty plasma implies that currents must close. 
Although the magnetospheric plasma maintains as equipotentials the magnetic field lines 
linking the ring plane to the ionosphere, it cannot carry currents across the field lines. Thus 
currents that flow through the base of the plasma cloud must close in Saturn's ionosphere. 
In the wake of the cloud, the levitated dust rapidly combines with the positive ions that 
balance its charge, leaving the spoke trail behind the cloud non-conducting. 

2.1. The model 

We analyse a simple model that captures the relevant features of the system of a dusty 
plasma cloud, ionosphere, and magnetosphere. A 2D strip of material in the xy-plane at 
z = represents the dust grains (plus neutralizing ions) at the base of the plasma cloud 
at a given time. 1 Infinite sheets of a different material at z — ±L take the place of the 
ionosphere. The magnetosphere consists of a vertical magnetic field of strength B embedded 
in a massless plasma that maintains the field lines as equipotentials. All calculations are 
done in the rest frame of the strip relative to which the ionosphere moves with velocity 

1 Defined in this way, the strip resembles a strip of metal, in which the dust grains are 
analogous to the ion lattice and the neutralizing ions are like the electrons which balance 
the charge on the lattice. 
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v = (0,t>, 0) (see Fig. 2). Thus x represents the radial direction in the rings, and y is 
azimuthal. 

We wish to determine the drift velocity of the plasma in and above the strip. This is 
equivalent to finding the horizontal electric field in the system since the drift velocity v p 
measured in the same frame as E satisfies 



2.2. Finding the drift velocity 

In the rest frame of a conducting sheet, the height-integrated current density J 
perpendicular to the magnetic field is given by 

J = ( xE + £^, (2) 

where a and £ are, respectively, the height-integrated direct (Pedersen) and Hall 
conductivities. Because the ionosphere consists of two sheets in parallel, its effective 
conductivity is twice that of a single sheet. 

Subscripts s and i are used to denote properties of the strip and ionosphere. The 
electric field in the rest frame of the ionosphere is related to that in the rest frame of the 
strip by 

Bxi B xs ~t~ 5 By S . (3) 

We assume, as GM implicitly did, that the currents are small enough so as not to 
significantly perturb the externally imposed magnetic field B, i.e. AtcJ/c <C B. We also 
take w<c. 

Then, using Eq. (2) for the currents in the rest frame of each conductor, and expressing 
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current conservation by 



J s + Ji = , 



(4) 



we solve the resulting simultaneous equations to give 



E, r . 
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(6) 



where E t = E s + Ej and cr t = a s + cr. 



The velocity v p at which the plasma cloud moves then follows from Eq. (1): 



EjEf + (TjCTt 
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(7) 
(8) 



where we have introduced b — B/\B\ — sign(S). 



The main message from Eqs. (7) and (8) is that |i> px |j kpj/| < v - 2 This finding 
contradicts the result obtained by GM (i.e. |i> px | = 20 — 70 km s -1 ), because in the region 
in which spokes appear, v < 1 km s^ 1 . The reasons for this discrepancy are detailed in §3. 

Another interesting limit is one in which either or both of the direct and Hall ionospheric 
conductivities are much larger than both components of the strip's conductivity. In this 
case, |i> px | <C v and v py ~ v; i.e., the plasma cloud basically corotates with the ionosphere. 



Because Hall conductivities can be negative, we can envisage the contrived case 2j s ~ 
— Ej, in which we can in principle have |t> px | > v. However, this requires the unlikely 
cancellation of the two unrelated Hall conductivities to a high degree of accuracy. Moreover, 
\vpx\ > v requires Ej > o"j, which is not the case for Saturn's ionosphere (see section 4). 
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3. Where GM Erred 

GM neglected several terms when solving equations analogous to our Eqs. (l)-(4). 
The components of the height integrated currents in both the strip and ionosphere as a 
function of the electric fields in their respective rest frames satisfy: 



J xs = [a s E xs ] - [bZ a E ya ] (9) 

Jy S = [<J a E ys ]+bll s E xs (10) 

Jxi = &iE X i — [bY>iEyi] (11) 

Jyi = &iEyi + [bT,iE x i] . (12) 



Terms left out in GM's analysis are enclosed in square brackets. Since the strip has 
negligible direct conductivity, neglecting terms proportional to a s does no harm. The Hall 
current in the ionosphere should be included since the Hall conductivity is comparable to 
the direct conductivity in the ionosphere. However, this neglect is not a large source of 
error. The serious omission is that of the strip's Hall current from J xs , as indicated in Eq. 
(9). 

The above equations, with the bracketed terms omitted, yield the incorrect result 

v px = v—b, (13) 

from which \v px \ > v follows for |E S | > <7j. 

Physically, GM's error is described by their incorrect statement that "the motion of 
the plasma cloud does not constitute a perpendicular (to B) current as the plasma cloud 
is charge neutral" . Because the dust grains are negatively charged the plasma must have 
a net positive charge. The azimuthal electric field set up by the differential motion of 
the plasma and the dust grains does cause the plasma to drift radially, but because the 
plasma is net positive, this constitutes a radial Hall current. This radial current closes in 
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the ionosphere, and thus modifies the radial electric field in both the ionosphere and the 
strip. But the radial electric field is responsible for driving the azimuthal Hall current in 
the strip, and this modification is not taken into account in GM. They therefore incorrectly 
fix v py = v, leading to a severe overestimate of v px . All these factors are accounted for in 
the self-consistent calculation outlined in §2. 



4. Recalculation of Drift Velocity 

We recalculate the drift velocity of a plasma cloud from Eqs. (7) and (8). 

The direct conductivity in the strip is very small because the dust grains have low 
mobility and the plasma is tied to the magnetic field lines. The Hall conductivity, which 
results from the E x B drift of the positively charged plasma, is correspondingly high, of 
order 

(14) 

where N p e is the height integrated charge density of the plasma. 3 Using the approach of 
GM, we estimate the height integrated charge density to be of order 10 1 esu cm -2 , which 
gives S s ~ — 10 14 cm s" 1 . 4 

Saturn's ionospheric conductivities vary with latitude and with time, with typical 
dayside values of the height-integrated direct conductivity being a { ~ 10 12 — 10 13 cm s _1 , 
and nightside values about 100 times smaller (Cheng and Waite 1988). The Hall 
conductivity is about an order of magnitude smaller, e.g. for the auroral region we have 



3 The Hall conductivity is negative because it is defined as positive for the commonly 
encountered case in which electrons are the dominant current carriers. 

4 This number may be spuriously high, because the approach of GM gives more charge on 
the dust grains than was originally present in the plasma. 
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<7j = 5 x 10 13 cm s -1 and Sj = 8 x 10 12 cm s _1 (Atreya al. 1983). The direct current 
is carried predominantly by protons, and the Hall current by electrons, which suffer 
fewer collisons per gyroperiod than the protons, so E x B-drift more freely. Collision 
frequencies of electrons and protons in Saturn's ionosphere are smaller than the respective 
gyrofrequencies, so the Hall conductivity is smaller than the direct conductivity. 

We substitute into Eqs. (7) and (8) the typical dayside values 2<jj = 1 x 10 13 cm s" 1 and 
2Sj = 1 x 10 12 cm s _1 (the factors of two account for both hemispheres of the ionosphere, 
although we note that we do not expect north-south symmetry of conductivities). We then 
obtain 

v px = 6 x l(rV v py = -2 x 10~ 3 v, (15) 

where we have used b = —1 as is appropriate for Saturn, and where v < 1 km/s. With 
these velocities, spokes will not form quickly or radially. The strip Hall conductivity is high 
enough to drag the plasma column in an almost Keplerian orbit. 



5. Discussion 

In response to this paper, Morfill & Thomas (2005) revisit the GM model for spoke 
formation. They provide further details about dust charging and plasma cloud structure. 
These details are however irrelevant to the criticisms raised here, and make no difference 
to the fact that, just as the radial electric field produces an azimuthal Hall current, an 
azimuthal electric field will produce a radial Hall current. These currents and electric fields 
must be treated self-consistently, resulting in the conclusions reached in this paper. 



6. Conclusions 



We have studied the physical situation in which a strip of conducting material moves 
between two parallel sheets of a different material, to which it is joined by perpendicular 
magnetic field lines (Fig. 2). The plasma on these field lines will drift in the plane of the 
strip, both parallel and perpendicular to the relative velocity vector. We have shown that 
the magnitude of this drift velocity cannot exceed that of the relative strip-sheet velocity. 

Application of this limit to the most popular model for spoke formation demonstrates 
that the model rests upon a gross overestimate of the velocity at which a plasma cloud can 
drift. 
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Plasma moves on: 
levitates more dust 
in new location 
to maintain E field 
and drift 



Extent of dust+plasma cloud 




Dust negatively charged: 
plasma has net positive charge 



Current closes at back and front of cloud 



ExB drift of all plasma particles: 
gives net current 



Fig. 1. — The dusty plasma cloud, viewed from above in the local rest frame of the ring 
particles. 




Fig. 2. — The ionosphere-strip configuration described in §2.1 



